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ABSTRACT 

We consider the formation of satellites around the Pluto-Charon binary. An 
early collision between the two partners likely produced the binary and a narrow 
ring of debris, out of which arose the moons Styx, Nix, Kerberos and Hydra. 
How the satellites emerged from the compact ring is uncertain. Here we show 
that a particle ring spreads from physical collisions and collective gravitational 
scattering, similar to migration. Around a binary, these processes take place in 
the reference frames of ‘most circular’ orbits, akin to circular ones in a Keplerian 
potential. Ring particles damp to these orbits and avoid destructive collisions. 
Damping and diffusion also help particles survive dynamical instabilities driven 
by resonances with the binary. In some situations, particles become trapped 
near resonances that sweep outward with the tidal evolution of the Pluto-Charon 
binary. With simple models and numerical experiments, we show how the Pluto- 
Charon impact ring may have expanded into a broad disk, out of which grew 
the circumbinary moons. In some scenarios, the ring can spread well beyond the 
orbit of Hydra, the most distant moon, to form a handful of smaller satellites. If 
these small moons exist. New Horizons will hnd them. 

Subject headings: Kuiper belt: general - planets and satellites: formation - 
planets and satellites: rings - planet-disk interactions 
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Introduction 


With its rich system of satellites, the Pluto-Charon binary seems a miracle of plan¬ 
etary dynamics. Separ ated by about 17 Pluto radii (1-Rp ~ 1200 km; lYoung fc Binzel 


1994 lYoung et al.l 120071). Pluto and Charon have masses of Mp 1.3 x 10^° g and Mq 


1.5 X 10^"^ g fjChristv fc Harringtonlll978l: iBuie et al.l 120061) . Their moons are not far away: 


Styx, Nix, Kerberos and Hydra are all packed between 37 Rp and 60 Rp, on coplanar. 


(Weaver et ah 2006: S 

howall 

ter et al. 

2011 

2012. 

2013; 

Buie et al. 

2013: 

Brozovic et al. 

2015; 

Showalter & Hamiltoi 

] 2015 

). The total mass of the satellites is a small fraction (< 0 

.02%) 

of the binary’s mass, roughly 3 x 10^® g to 3 x 10^^ g 

Buie et al 

2006: 

Brucker et al. 

2009; 

Showalter et ah 

2011, 

2012; 

Youdin et al. 

2012; 

Brozovic et al. 2 

0T5h. 

Resonances with the 

e satellite svstem on 

binary and the gravitationa 
the verge of chaotic disinteg 

. effects of the moons on each other keep th 

ration (j" 

rholen et al 

2008 

Siili & Zsigmonc 

2009: Winter et al. 

2010; 

Peale et al 

20i; 

.: Youdin et al 

.2012 

: Showalter & Hamilto 

n 2015 

). How the satellites 

Tallenge for theories 

arrived at their fragile orbits around Pint 

This delicate orbital architecture mak( 

3 and Charon is a mystery. 

3S the Pluto-Charon svstem a 

of circumbinarv planet formation fe 


doriwaki & Nakagawa 

20041: ( 

Duintana & Lissauer 


20061: IPierens fc NelsonI 120071: lRafikovll2013l: iBromlev fc Kenvon 


201511 . The most promising 


scenario for the origin of the binary itself is that two dwarf pl anets had a grazing colhsion 


1992 


leaving thern intact but bound at a small separation of 4-5 Rp f McKinn^ 1989 : Stern 
Canupl 120051. 1201111 . Pluto and Charon then tidally evolved, moving out to their present 


positions, locked synchronously together with a 6-day orb ital period fiFarinella et ahl 11979 
Dobrovolskis et ahl 119971: iPealel Il999l: ICheng et ahl l2014bl) . The impact likely produced a 


dynamically hot ring, confined to a region no more than a few times the binary’s initia l 
separation, inside the present-day orbit of Styx, the innermost satellite (jCanupll2005l. 1201111 . 


2006; 

Canup 

2005 

). After settling at orbital distances of 5-30 Pluto radii ( 

Canup 

2005). 


the debris circularizes into a narrow circumbinary ring. Over time the ring spreads radially 
outward, and small debris particles within it grow through coagulation. The result is a 
satellite system that has e xpanded outward by a factor of three or more in orbital distance 
( Kenvon Sz BromlevI 2014 ). A major uncertainty in this picture is how the ring — either 
with fully formed satellites or smaller debris — can spread to the moons’ current positions. 
A compelling idea is that the moons formed quickly a nd migrated outward, trapped in 


resonances that expanded as the binary tidally evolved (jWard fc Canupl 1200611 . However, 


this mechanism depends on details o f the tidal evolution, and no single model can place all 


of the moons in their present orbits flLithwick fc Wull2008l: ICheng et al.ll2014all . 
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Despite difficulties with the lWard fc Canupl (120061) resonant migration scenario, it has an 
important implication. Even if satellites were placed by hand in their present orbits just after 
the giant impact, then at least sor ne would be lost when unstable resonances swept outward 
along with the expanding binary ( Lithwick fc Wu 2008 : Cheng et ah 2014alf5) . Collisional 
damping and diffusion might mitigate the effects of the resonances, but only if satellites 
evolv e in partnership with a reservoir of smaller particles to damp them (jWalsh fc Levison 
20151 ). Therefore, the key to putting the moons into their current orbits is likely coordination 
between growth of the satellites, the depletion of the smaller debris, the spreading of the ring, 
and the tidal evolution of the binary. The primary motivation for this work is to understand 
this interplay. 

Here, with the Pluto-Charon system in mind, we investigate circumbinary ring dynam¬ 
ics to see how a compact ring of growing particles evolves. Thro ugh analytic estimates 
and numerical calculations, w e examine th e role of particle viscosity (ICook fc Franklinlll964J: 
Goldreich fc Tremaind Il978[) and gravity flLin fc Papaloizoul Il979al: IWardI 119971) . We high¬ 
light the dependence of the spreading rate on the state of the ring particles, particularly 
their random velocities wit hin the ring, and we consider the effects of in situ satellite growth 
(iKenvon fc BromlevI 120141) . In developing simple parameterized models, we offer pathways 
from an initially compact ring around a short-period binary to the more extended system 
of moons we see today. Some models predict the formation of more small moons at orbital 
distances beyond Hydra. 

In this paper, we outline the main physical phenomena that are important to the for¬ 
mation of a circumbinary satellite system like Pluto-Charon’s (§2). Then we give a general 
introduction to ring dynamics (§3), including the effects of a central binary on the ring 
(§4). Turning our analysis to Pluto-Charon, we describe models of ring evolution (§5), 
and conclude with sor ne predictio ns for New Horizons (§6), now rapidly approaching the 
Pluto-Charon system fjSternl 120081) . 


2. Overview and Context 

The formation of the Pluto-Charon satellite system is a challenging problem involving 
a range of physical phenomena. The grazing collision scenario for the formation of the 
binary yields only a narrow, compact ring that somehow spreads, either from interactions 
between ring particles or with the binary. The time scale for these processes is uncertain. 
In some circumstances, the ring spreads rapidly through gravitational stirring and physical 
collisions. In other situations, the ring spreads slowly as a result of collisional damping 
with only long-range, weak gravitational interactions to drive the expansion. In both cases. 
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resonance excitations from the binary sweep throngh the system and potentially destabilize 
satellite orbits. One of onr primary goals is to understand how the ring expands as the orbit 
of the binary evolves. Figure d] illustrates the main issues, along with the overall layout of 
the Pluto-Charon system. The following list provides more detail about characteristics of 
the binary, the ring around it, and the physical processes that guide the emergence of the 
satellite system. 


The circumbinary environment: the inner cavity. The Pluto-Charon binary 
provides a time-va rying potential that complicates the nature of orbiting satellites 


(ILee fc Pealell2006l) . Circumbinary orbits near Pluto and Charon are largely unstable. 


of Styx (e.g., the light gray shaded region in FigurefH see also 

dolman & Wiegerl 

1999; 

Musielak et al. 2005: Pichardo et al. 

2005 

Doolin & Blundell 

2011) 

. Carved out by in- 

stabilities at overlapping resonances 

(Wisdom 19801) and bv resonant-driven tora 

ue ex- 

change (e.g.. Goldreich & Tremaine 

1980; 

Mever-Vernet & Sicardv 

1987: Ward 

1997), 


this cavity sets the inner boundary of the circumbinary ring. Due to uncertainties in 
this complex environment, we consider two approaches to the boundary: (i) loss of ring 
material by ejection and (ii) a balance between viscous inflow and the outward torque 
from Charon. 

Resonances. Aside from defining the edge of the inner cavity, unstable resonances lie 
withi n the inner cavity and (for eccentric binarie s) at isolated larger orbital distances 
(e.g., IWard fc Canupl 120061: ICheng et al.l l2014bl) . Current analyses favor a circular 


Pluto-Charon orbit; resonances beyond the cavity’s edge at the 3:1 commensurability 
are sta ble (e.g., Figu re [3l below). Initially, the binary probably had significant eccen¬ 
tricity (jCanupl 120051) : resonances at higher order commensurabilities, include those at 
4:1, 5:1, and 6:1, were then unstable. Here we explore the effect of these resonances 
in the context of ring dynamics f 1 14.II) and show that sufficiently strong collisional 
damping stabilizes particles within these resonances. 

Stable circumbinary orbits. Despite the resonances, the binary has little impact on 
muc h of the orbital domain outside the in ner cavity. As demonstrate d by analytic the¬ 
ory (Lee fc Pealell2006l:Ifjenng fc T;eell2013f) and numerical experiment (IBromlev fc Kenvon 
20151) . non-resonant orbits experience no secular excitation or long-term torque ex¬ 
change with the binary (§!]). Particles tend to follow trajectories that ebb and flow 
with the orbital motion of the binary. When swarms of particles experience orbital 
damping, they settle on a family of “most circul ar” paths that do not in tersect, just 
like circular paths around a single central mass flYondin &: KenvonI 120131) . For a ref¬ 
erence frame tied to a most circular orbit, the local dynamics within a ring around a 
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binary is very similar to dynamics around a single central mass. After describing these 
features O, we take advantage of them when modeling ring evolution (^. 

Ring stability considerations: total mass and particle size. We also consider 
the gravitational stability of a ring with velocity dispersion n, orbital frequency fl, and 
surface density S orbiting Pluto-Charon. When ring particles achieve a balance be¬ 
tween gravitational stirring and collisional damping, the velocity dispersion is roughly 
the escape velocity of the largest particles in the ring. Dehning G as the gravitational 
constant, the Toomre sta bility criterion, > tt GS, establishes stable conhgurations 
of ring parameters (e.g., IChiang fc YoudinI l2010h . Adopting a ring radius of 20 Rp, 
width of 5 Rp and mass, 3 x 10^° g, comparable to the mass of the known satellites, 
the ring is stable if the velocity dispersion exceeds the escape velocity of a 10-meter icy 
particle, roughly 1 cm/s. Unless mergers produce larger particles which gravitationally 
stir the smaller particles up to their escape velocity, a ring composed of smaller parti¬ 
cles is unstable. In most of the models discussed here, we assume a sea of particles with 
radii of 1 km; rings composed of these particles are stable up to masses of 4 x 10^^ g, 
roughly two orders of magnitude more massive than the binary’s satellites. 

Ring dynamics. We envision the spreading of the ring from its early compact conhgu- 
ration through particle interactions. When particle velocities exceed their escape speed, 
either from stirring or because th ey are in the Roche zone of the central mass, physica l 
collisions dominate the dynamics (ICook fc Franklinlll 964J: iGoldreich fc Tremainelll9781 1 . 
Then, an effective collision-driven viscosity drives the evolution of the ring. Outside 
the Roche zone (which is only about 3 Rp for Pluto) collisional damp ing reduces the 


relative particle speeds uii til gravitational interactions are important flHornung et ah 


19851: IShu fc Steward Il985[) . Here, we provide a theoretical prescription for how these 
longer range interactions cause the ring to evolve. As illustrated in Figure [H collision- 
ally driven viscosity, with stirring by larger embedded particles, leads to fast migration. 
Purely gravitational spreading is comparatively slow. 

Evolution of the binary. Tidal evolution of the binary lies at the heart of all these 
issues. After the grazing impact, the binary has a tighter, more eccentric orbit. Tides 


circula rize and expand the system into the current circular orbit (e.g., ICheng et ah 
2014aJ). Our analysis considers how tidal evolution drives a set of sweeping mean 


motion resonances ( e.g.. Figure d]) which lead to the dynamical ejection of satellites 
i Cheng et ah 2014bh . By including collisional damping and ring viscosity, we begin to 
show how growing satellites or their precursors experience resonant trapping, allowing 
them to migrate outward with the binary’s tidal expansion. Alternatively, if the ring 
spreads quickly, we consider how satellites might survive, or at least reform, after the 
resonances sweep past. 
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While the circumbinary nature of Pluto-Charon satellites make them unique in the So¬ 
lar System, key physical processes governing their formation are also important to other 


cously spreading circumprimary disk (see Lunine & Stevenson 

1982|: 

JanuD & Ward 

2 OO 2 I: 

MosQueira & Estrada 2003a b: Mosaueira et al.ll2010l: Sasaki et a 

. 2010 

Ward & Ganur 

J 2 OIO: 

Ogihara & Ida 2012 

). Although gas dynamics is crucial to this evolution, our analysis ad- 


dresses aspects of the spreading of large particles uncoupled from the gas. With its gas 
free environment, the Pluto-Charon sate l lite system is more akin to Saturn’s rings (see 


Canup fc Espositolll995l: iPorco et al.ll2007l:ICharnoz et al.ll2010l) . Despite clear differences in 
the location of small particles relative to the Roche limit (e.g., mostly inside for Saturn and 
entirely outside for Pluto-Charon), both systems are driven by pure particle dynamics. To 
extend the theory of Saturn’s rings to a ring of solids orbiting Pluto-Charon, we focus next 
on the physics of particle rings. 


3. Ring dynamics 


Planetary rings provide an excellent theoretical proving ground for planet formation. 
Saturn’s rings famously became a subject for James Clerk Maxwe ll, and inspired a pow¬ 
erful theory for evolution w ithin the Roche limit of the planet flCook fc FranklinI ll964J: 
Goldreich &: Tremainel Il978l) . Inside this radial distance from Saturn, tidal forces inhibit 
particle growth. Thus, the theory concentrates on the radial spreading and velocity evolu¬ 
tion of small, indestructible particles with negligible gravitational interactions. Beyond the 
Roche limit, rapid growth into small moons is a likely outcome; radial sprea ding and grav¬ 
itational scattering within the rings remain secondary fjCharnoz et al.ll2010l) . For a ring of 
particles orbiting Pluto-Charon outside the Roche limit, gravitational interactions between 
large particles add to the spreading from small particles. Furthermore, dynamical spreading 
and particle growth occur simultaneously. To understand the evolution of this material, we 
require a more general treatment. 

Here we establish a theory for describing the ma jor influences on the evolution of a 
planetary ring. Following iGoldreich fc Tremainel (1197811 . we discuss the role of physical col¬ 
lisions under the assumption that the particles rebound off one another as they collide, 
changing their trajectories but n ot their masses or radii f ^3.2p . Inelastic collisions allow 
damping of random speeds (e.g., (Bridges et al.lll984l) and convert orbital en ergy into ran¬ 
dom motions, yielding an effective viscosity that s preads material radially fiJeffrevsl Il947l: 


Gook fc Franklin 1964: Goldreich fc Tremaine 1978 


below). We also consider the role 
of mutual gravitational interactions. Gravitational scattering affects the equilibrium veloc- 
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ities that enter in the viscosity flHornung et al.lll985l: IShu fc Steward Il985l: iGoldreich et ah 


20041) . a s well as the radial spreading of particles from collect i ve interactio ns similar to mi¬ 


gration fiLin fc Papaloizoulll979at IGoldreich fc TremaindllOSOl: IWardlll997l) . Finally, we give 
a prescription for the evolution of a ring in which particles grow. By incorporating all of 
these phenomena together, we construct simple models for ring evolution and lay the foun¬ 


dation for including this physics into full evolu tion codes fe.g.. IWeidenschilling et al.l 119971: 


Kenvon & Luu 

1998; 

Kenvon & Bromlev 

2014) 


3.1. Background 


Disks or rings around a central mass are a natural outcome in many astrophysical 
systems. Gonservation of angular momentum, along with dissipative processes (e.g., radiative 
cooling of gas or in elastic collisions between particles) drive material towar d these flattened 
conhgurations fe.g.. lLvnden-Bell fc Pringlelll974 iBrahidll 9761: |Pringlelll98lh . In this section, 
we examine a ring of particles around a central point pass in the limit where the self-gravity 
of the ring is not important. We begin by establishing some basic physical scales and 
relationships between physical properties of the ring. 

To model the ring, we hrst assume that it is made of identical particles with hxed radius 
r, mass m, and density p. A particle at an orbital distance a has an orbital period 


= 27rVaVGM, 


( 1 ) 


and Keplerian velocity 


= \/GM/a, 


( 2 ) 


where M is the central mass and G is the gravitational constant. We assume that the orbits 
of the ring particles are approximately circular and have instantaneous speeds close to Vj^. 

For orbits near the central mass, tidal forces are strong. The Roche limit gives a critical 
value. 




l7j 


(3) 


inside of which m assive particles are unable to hol d themselves together by their bulk strength 
and self-gravity I Holsapple fc Michel 2006 . 2008). 


When particles are wel 


outside of the Roche limit, they interact gravitationally with 


other neighboring satellites (ISalmon et ahll2010t iGharnoz et al.lDOllt iRosenblatt fc Gharnoz 



































































2012|). The gravitational range is characterized by the Hill radius, 


r^= a 


m \ 1/3 


/ m \ 


( 4 ) 


The particle’s Hill velocity, 


Ujj = a^/Gm/r^, 


(5) 


is the speed of a low-mass neighbor orbiting the particle at the Hill radius. Neighbors passing 
by at much faster speeds are undisturbed; slower objects are stirred to at least this speed 


see 


Goldreich et al. 


2nn4 . 


The fastest speed to which a p article can stir a less massive neighbor is approximately 
the escape speed at its surface (see iGoldreich et al.ll2004l) : 


= \/2Gm/r = ^/SnGpr'^/S. 


( 6 ) 


During an encounter, if the relative speed of an object and the satellite is small compared 
to Vesc, gravity influences the outcome. When the relative speed exceeds Vesc, gravitational 
interactions are much less important than physical collisions. 

In a swarm of identical particles, the escape velocity provides an estimate of the char¬ 
acteristic speed of particles, v, as measured in a local Keplerian frame that follows a circular 
orbit in the ring midplane. The actual value of v, which we associate with “random mo¬ 
tion” in this reference frame, depends on collision outcomes and physical conditions in the 
disk, such as the number density of particles n(a), the surface (mass) density S(a), and the 
vertical scale height h{a). These properties of the disk are inter-related. 


n(a) ~ 
h{a) ~ 


E 27rE 

_ _ 

2 hm mvT^ 
dvr ’ 


(7) 

( 8 ) 


where we assume a vertical speed v^, ~ v/2 (e.g., lOhtsukil Il992l. and references therein). 
Thus, the microscopic particle velocities and the broader structure of the disk are entwined 
throughout its evolution. 


3.2. Particle collisions 

Physical collisions between particles hold the key to the evolution of systems of satel¬ 
lites, moons, and planets. Gollisions not only result in fragmentation or mergers; they also 
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damp random particle motions and reduce collision speeds. In the balance lies the differ¬ 
ence between growth or destruction of satellites in a ring. When particles collide, their 
material properties determine the outcome. For icy bod ies, shattering or fragmentation oc¬ 
curs at speeds of cm/s flKenvon &: Bromlevll2004 and references therein). At slower 

speeds, bouncing and sticking can occur. Because high-speed collisions produce smaller bod- 
ies that can dynainically cool larger ones, collisions are essential to mergers and growth (e.g., 
Youdin fc Kenvon 2013 ). 


A first step toward assessing the role of collisions is to estimate T^oi, the characteristic 
time between collisions for a single particle. If particle velocities are comparable to or larger 
than their escape velocities, simple kinetic theory (“ring ”) gives fe.g.. ISafronovlll969l: iLissauer 


19871: IWetherill &: Stewartlll993 


Goldreich et al.l 120041) 




col 


pr rj. 

2x/27rE 


(9) 


n\/2v'KAr‘^ 

where the cross-sectional area of 7r(2r)^ is appropriate to hard spheres, and the pairwise 
velocity \/2n is eliminated in the ri ghtmost expre ssion because of its relationship to the 
scale height, h (Equation ([H])); see lOhtsukil Il992l ). When the relative speeds fall below 
the typical escape speed, gravitational focusing i ncreases the effective cross-section and 
reduces the collision time ( Wetherill fc Cox 198,4 Wetherill fc Stewart 1989 : Spante et ah 
1991 : Goldreich et ah 2004 : Youdin fc Kenvon 2013, see 1 13.4.21 below). 


If collisions are inelastic, relative particle velocities are damped. Absent other effects, 
we can identify a damping time scale in terms of the evolution of random kinetic energy: 


dt 


1 


■ damp 


Thus, the damping time is 


Z 


■ col 


damp 


/lo 


( 10 ) 


( 11 ) 


where /loss is the average fractional loss in kinetic energy per collision. Since the random 
speed V is measured in the local Keplerian frame, the damping time characterizes how colli¬ 
sions circularize the orbits of ring particles and flatten the ring. 

The paramet er /loss carries the details of collisions between ring particles. Following 
Porco et al.l (120081 ). we adopt simple parameterizations of the coefficients of restitution, which 
set the ratio of the particle speeds before and after a collision. Treating the normal (N) and 
transverse (T) directions separately, these coefficients are 


Cat — 

ex = 


mm 


{vn/v.) ,1 


( 12 ) 

(13) 


0.9, 
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where vn is the normal component of the relati ve velocity. Experi r nents snggest = 
0.01 cm/s and b = —0.14 for spheres of water ice flSupnlver et ahlll995l: [Bridges et ahlll984 


Porco et ahN2008h . Averaged over all impact parameters, the loss parameter for a given 


collision speed u is 


flo 


[U 


= 1 - 


{e% + e| 


)/ 2 , 


(14) 


if collisions are random as in a 3D ideal gas. Averaging over the distribntion of relative speeds 
yields /loss (Eqnation (ITTIB . In a simple model with a Maxwell-Boltzmann distribntion of 
relative speeds, 


_ u^dufiossju) exp(-3nV4n^) 
u^du exp(—3-u2/4n2) 

/ \ - 0.28 

^ 0.60 — 0.44 (—j [6 = —0.14,n 3> n*]. 


(15) 

(16) 


In this case, the energy loss for n S> n* is limited by grazing collisions and the valne of er- 
Unfortnnately, this valne is not well constrained by experiment. 


3.3. Ring viscosity 


While physical collisions sap random kinetic energy from the ring, they also add to 
it. By redirecting their trajectories after rebou nding, collisions convert orbital motion into 
random velocities (IGoldreich fc Tremaine! 1 19781 ). Thus, the eccentricities and inclinations of 
ring particles grow; the entire ring broadens in semimajor axis. The diffusive spreading of 
ring particles can be expressed formally in terms of a viscosity. 


Dad 


IjT 
0.46- 


r 


2ti 1 + ' 

where the optical depth of particles in the ring is 

3S 
4pr 


r 


m 


(17) 


(18) 


fjCook fc Fra.nk]inlll964l: IGoldreich fc Tremainelll978l. Il982[) . To illustrate how the viscosity 
i^rad is associated with diffusive spreading of the ring, we consider a particle as it rebounds 
during a collision. If we are oblivious of the details, the collision outcome appears stochastic, 
yielding a random step with a radial component that is directed either inward or outward 
with equal probability. The magnitude of this step is roughly the particle’s epicyclic radius. 


|(5a| ~ a- 


( 19 ) 
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Viewing this outcome as part of a random walk, with steps taken at the collision rate of 
1/T’coi (Equation ([9])), we can infer an effective diffusion coefficient, 


^rad 


|(5a| 










col 


v^Tcoi (27ra/Tj^)2 • (prT^/Z) Anpr 


37r ’ 


( 20 ) 


where the rightmost expression is the low optical depth limit of Equation flTTl) to within a 
factor of ~ 2. After n steps (collisions), the typical radial displacement grows to 


\5a 


(Vi 


\fiA ~ a—^t/Tcoh 


( 21 ) 


with time t set to nTco\- We can take advantage of this formula in simulating particle 
diffusion. Using a random number generator we can send particles on randorn walks in 
orbital distance to mimic the behavior of a viscous ring flBromlev fc Kenvonll2013l. also ( 14.ll 
below). 

Over many collisions involving an ensemble of particles, the ring broadens by some 
radial distance |Aa|. The time scale for this spreading is 


Z 


spread 


|Aa^| 

Uad 


( 22 ) 


This relationship explains the sl ow spreading of Satu rn’s rings. Doubling the annular width 
of the A ring requires ~ 1 Gyr fjCharnoz et al.ll2009h . Slow spreading is a major challenge 
for models in which an impact ring ar ound the Pluto-Charon bi nary expands in time to form 
the moons at their present location flKenyon fc Bromley! 120141 ). We return to this issue in 


Because collisions move particles to regions with different Keplerian velocity, they tend 
to pump up the relative sp eed of the particles. If collisio ns are completely elastic, the relative 
velocities grow unchecked flGoldreich fc Tremainelll978h . For example, a particle in a sea of 
objects on circular orbits might scatter slightly inward or outward after a collision, so that it 
has a new speed relative to its new neighbors, |Ar>| ~ e|Aa|. Subsequent collisions cause that 
particle to walk randomly in velocity, as well as in radial distance. By the interaction 
the typical speed has grown to 




n) 1 2 


c\Av 


(n-l)|2 


~ C‘ 


n—1 


|An| 


(23) 


where c is a constant greater than unity that depends on the details of the int eractions (e.g., _ 

parti cle sizes, orbital distances). Formally, random speeds quickly diverge (e.g.. lGoldreich fc Tremaine 

1978h . 
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In reality, collisions involve some loss of kinetic energy, providing a means to balance 
this process. Coefficients of restitution that decrease with speed — and therefore cause more 
kinetic energy loss at higher speeds — allow a ring to reach an equilibrium. If the collisional 
damping is weak, rings equilibrate at high speeds; highly inelastic collisions lead to small 
random motions and nearly circular orbits. The details of this balance between dynamical 
heating and collisional damping are buried in the magnitude of the equilibrium velocity, v 
in Equation (1T7|) . 


3.4. Gravitational scattering 


In addition to physical collisions, pairwise gravita tional encounters can also cause ran¬ 


dom deflections and thus contribute to ring viscosity flHornung et al.l llQSSt iBarge fc Pellat 
1990l: IOhtsukilll999l: lOhtsuki et al.ll2002h . From Equation (120|) . we can surm ise how gravita- 


tional interactions affect Vr 


Gravitationa l stirr i ng and dynamical friction (jWeidenschilling 


19891: kenvon fc Luul Il998l: IStewart fc Idal 120001: iKenvon fc BromlevI l200lh drive the local 
velocity evolution; gravitational focusing makes the effective cross-section larger and reduces 
the collision time. Thus, stirring and focusing tend to increase the viscosity and drive the 
spreading of the ring. However, for rings in a steady state, stirring pumps up velocities until 
they exceed the escape speed of particles. Then, gravitational focusing becomes unimportant 
(see below) as collisions take over to regulate the velocities. Therefore, gravity’s main role 
in the ring viscosity is to help to maintain the equilibrium velocity established by collisions. 
After all, only three things matter to the spreading of the ring: how many particles are in 
the ring (S), how big they are (r) and how fast they are moving relative to each other (v). 

Here we review some general aspects of gravitational interactions between ring particles. 
For particles outside the Roche limit, separat ed by distances of a few H ill radii, we describe 
pairwise interactions as Rutherford scattering flLin &: Papaloizoulll979bl) . Working in a pair’s 
center-of-mass frame, we relate impact parameters to scattering angles. Translating these 
scattering outcomes to changes in orbital elements, we then derive how gravity between ring 
particles affects the overall evolution of the ring. 


3.4-1- Rutherford scattering 


Pairwise gravitational interactions between massi ve particles are described by Ruther¬ 
ford scattering. Following the traditional analysis (e.g.. lLin fc Papaloizonlll979bl l we consider 
two point particles with masses rui and m2 on approaching trajectories in their center-of- 
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mass frame with impact parameter b and closing speed v. Using the total and rednced 
masses 

m = rrii + m 2 and 6m = -, (24) 

mi + m2 

we write eqnations of motion for the pair as a one-body problem in terms of the pair’s relative 
position and velocity. We are typically interested only in the ontcome of these enconnters, 
described in terms of 

^ Gm , , 

0 = 2 arctan , (25) 

vm 

which is the angle between the initial and hnal relative velocities. Another qnantity that is 
important to scattering ontcomes is the distance of closest approach, derived from 


^2 _ ^2 _ ^ Gmr ^,_ 




(26) 


where Wesc is the escape velocity of an individnal particle of radius r. Between the distance 
of closest approach and the scattering angle, we distinguish between different outcomes: 
physical collisions (r is less than the sum of particle radii) are one possibility as described 
above. For purely gravitational encounters, we also have strong/random scattering (0 ~ 1, 
Uiose weak scattering (0 -C 1, h). We consider these regimes next, starting 

with strong interactions. 


3.4-2. Gravitational focusing and random scattering 


When gravitational encounters are close and deflections are strong, physical collision 
rates can be boosted and relative velocities can be randomized. In a pairwise encounter, if 
the closing speed v is small compared to the escape velocity of the largest particle in a pair, 
then the distance of closest approach r is signihcantly smaller than the impact parameter 
b (Equation ([26]). This effect is gravitational focusing. If we compare the cross section for 
physical collisions between identical particles to be, the maximum impact parameter that 
leads to a collision, then we get the focusing factor. 


7r(2r)2 


1 + 


2^2 

^esc 


(27) 


This factor quantihes the boost in the like lihood of physical collisions and in the collision 


rate (in the context of plan et formation, see lGreenzweig fc Lissauerlll990l: ISpaute et al.lll991 


Wetherill fc Stewartlll993l ). The overall rate of ring spreading is increased, with the viscosity 


parameter augmented by a factor /fo 
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Even in the absence of collisions, strong gravitational enconnters can affect the viscosity 
by randomizing velocities. When the impact parameter b is comparable to Gm/v'^, ronghly 
the “gravitational radins” of a particle, the scattering angle 0 covers the fnll range of pos¬ 
sibilities. Thus, if the impact parameter is stochastic, the scattered velocities are as well; 
close gravitational interactions act like elastic collisions. Then, to derive the spreading of a 
ring, we may adopt the formalism for collisions described in 1 13.31 Thus, in the optically thin 
limit, there is an effective viscosity that describes both physical collision s and gravitational 
scattering. (For rigorou s tre atments of diffusive s preading by gravity see lOhtsuki fc Tanaka 


20031. lOrmel et al.ll2012l and iGlaschke et al.ll20l4 i 


3.4-3. Collective gravitational effects: migration 


In addition to contributing to random motions, pairwise gravitational interact ions also 


act collectively to produce steady flows. Studies of planetary and satellite migration fiLin fc Papaloizou 
1979a : Goldreich fc Tremainel 19^: Ward 1997 1 show how swarms of particles contribute over 
relatively long distances, O(10)rjj to cause a steady radial drift of mass in the ring. 

To quantify this effect, we start with the Rutherford formula (Equation fl2^ ). and 
work only in the limit of small scattering angle 0 since we are interested in larger-scale 
collective effects, not small-scale random motions. As above, we consider two particles with 
masses mi and m 2 on approaching trajectories in their center-of-mass frame with impact 
parameter b and closing speed v. When the two particles are both orbiti ng a central body, 
scattering through an angle 0 tran slates to changes in orbital elements flLin fc Papaloizou 
1979bl: iGoldreich fc Tremainel 1198011 . To see how this process plays out, we first consider a 
pair of particles with extremely unequal masses. 

We suppose mi —)■ m and m 2 —)■ Sm m and that the more massive particle is on 
a circular path around a central object of mass M at an orbital distance a. Similarly, the 
lighter body has an orbital distance of a + q, with q/a -C 1, orbiting on a circle in the same 
sense. We can then write the denominator of equation as a product of the pair’s relative 
orbital speed, 

and the relative orbital angular momentum 


i ^ SmCtaq 


(1 -I- qjaf’C 


(29) 


where G is the orbital angular velocity of the massive body. The change of angular momen- 
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turn of the smaller object is 

= (1 — cos 0) £ -0^£, 


(30) 


in this small-angle, “weak scattering” limit. Combining these results and expanding in terms 
of g, the change in orbital angular momentum is 




16 G‘^m?5ma 
~9 


(31) 


This perturbation of the small satellite causes an equal and opposite change in angular 
momentum of the larger particles, Ai = —6i. The semimajor axis of the larger object also 
changes. For circular orbits, 

f ^ GMa. (32) 


Taking derivatives of both sides of this expression 
semimajor axis as 


6a ^ 


25i 

mVLa 


gives the change in the larger body’s 


(33) 


Applying Kepler’s third law {GM = 


6a 


32 m 6m aP 
9 g® 


9 g 

1 +-- + ... 
4 a 


[m S> 6m] 


(34) 


A similar analysis of scattering between equal mass objects gives 

256 m2 a® 


6a 


25 M2 g5 


, 27 g 

1 H-- -h ... 

20 a 


[equal-mass satellites] 


(35) 


The change in orbital distance per encounter per unit mass is 


6am{a,q) ^ 


8 Am a® 
M2 g5 


1 + B« 


(36) 


where A = 4/9 and B = 9/4 when m S> 6m, and A = 32/25 and B = 27/20 when the mass 
of the scatterers is equal. Thus, in general A = 0.4-1.3 and B = 1.3-2.3 (the factor of 8 in 
Equation (l36|) was chosen to keep these constants near unity). 

Equation fl36|) is valid only for a limited range of orbital configurations. The Hill radius 
of the larger satellite, which characterizes that particle’s domain of influence, provides a lower 
limit to the orbital separation g. Furthermore, for the weak theory to apply, an interaction 
with separation g must yield small orbital perturbations. This condition is violated if the 
orbital separation is within 

« 3.5r„ (37) 
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since the outcome of such encounters is that the pair chaotically cross orbit^ f Glad man 
19931) . The theory breaks down at large orbital separations, where the cur vature of the 


satellites’ orbits plays a role in the interactions flBromlev &: Kenvonl 1201 Ibh . Combining 
these two constraints gives 




(38) 


which is an approximate condition for Rutherford-style weak scattering. 


Migration theory uses these conclusions to predict the drift rate of a satellite in a particle 
disk by tracking the rate of encounters along with the orbital displacement per encounter 


(e.g., Ilda et al.l 120001) . This same strategy applies to a ring of equal-mass particles. The 
interaction rate between pairs of particles at orbital separation q is set by their synodic 
period: 


Tsyn(®) q') 


2alh 


3|g| [l-5g/(4a)] ’ 
Integration over the disk gives the total drift rate, 

S(a q)6am{a, q) 


(39) 


n 

^grav 


27r(a + q)dql-i{q)- 


Tsyn q^ 


(40) 


where l-i^q) is a high-pass filter, set to zero at small scales (e.g., |g| < 3rjj) and unity elsewhere. 
Thus the filter eliminates contributions to the drift rate on scales where the weak-scattering 
theory breaks down and where interactions are dominated by random motions. 


With the weak-scattering result for 5a, and a Taylor expansion in terms of the radial 
distance variable g, the drift rate becomes 


^grav 


24R7ra2E m a 
M MT 


sgn(g)'H(g) 


dq 
a q^ 



V 4 E da J 


(41) 


where the constants A and B are of order unity and depend on the mass ratio of the drifting 
particle and other objects in the disk (see Equations fl3T)l and ([35])). 

In a dynamically cold ring — where the scale height h is smaller than the Hill radius 
of the ring particles — the minimum orbital distance between particles for which steady mi¬ 
gration is possible is ~ g^j^^ (Equation (1371) b Inside this distance, orbits are chaotic and 
can cross, enabling physical collisions. Outside this scale, orbital evolution is diffusive. The 
ring of particles expands until particle separations become very large (e.g.. Equation (155]) b 


^For exactly circular orbits in the extreme mass-ratio case, weak deflections can cause orbits to cross 
when the orbital separation is less than 
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If the ring is dynamically hot, with random particle excursions comparable to the scale 
height h, then ^ h. By setting 'H(g) = 0 when 1^1 < and unity elsewhere, we obtain 
a general expression for the radial migration of a ring particle: The radial drift rate of a 
particle embedded in the ring is 


^grav 


247ra^S m a ^ f ^ 1 a d'L\ 

■^min / 


(42) 


where 

= max(g^j„^, h) or ^ (Ibr^ + h?). (43) 

The rightmost expression just serves to enforce the minimum value of q^^^ as a smooth func¬ 
tion of h or Tjj. In the expression for the drift rate itself, the quantity i? — 1/4 represents 
a systematic inflow of geometrical origin in a disk with uniform surface density, while the 
logarithmic derivative corresponds to a radial flux in a non-uniform disk analogous to diffu¬ 
sion in a viscous ring. Together these two terms dehne an equilibrium surface density with 
no radial flux (din S/d In a = 1/4 — 5). Near the inner edge of the ring, where the radial 
gradient of S is large, material flows inward. Similarly, the ring spreads outward near the 
outer edge. 


3.5. Ring evolution 


Astrophysical rings continuously evolve with time. Physical collisions and gravitational 
interactions produce diffusive spreading, growth, and destruct ion. Inside the Roche limit. 


simplifying assump tions, such as requiring hxed particle size flGoldreich &: Tremaine! 1 1978 


Salmon et al.N20101 1. allow elegant solutions to this evolution. Here, our goal is to track 
simultaneous growth and viscous spreading. We take advantage of the idea that pairwise 
interactions drive all the dynamics. In a steady ring, most interactions involve scattering and 
are neither destructive or fully inelastic (all collisions lead to perfect mergers). Then, the local 
dynamics (viscous or gravitational spreading) can be handled by approximating particles as 
having a hxed size over a certain short time scale to track the dynamical evolution of the 
ring. Particle growth can then enter as a separate step in the evolution, in staggered fashion, 
in the spirit of a st andard leap-f rog technique. This approach also allows us to use standard 
tools of the trade flPringlelllQSlI) to calculate the ring evolution. 
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3.5.1. Steady-state theory of ring evolution 

We describe the evolution of a ring in terms of its surface density. The e volution equation 
follows from conservation of mass and angular momentum flPringlel 1198111 . Integrated over 
the direction perpendicular to the ring plane, these conservation laws are 


dt 


d 


+ 


1^ 

a da 

1 d 


(aaS) = 0, 


avails) = 


1 dL 


(44) 


(45) 


dt '' ^ ' ada ^ ^ 2Tia da ’ 

where a is orbital distance, hi is orbital angular speed at distance a, and L represents the 
torque on an annulus of radius a and width Aa. These equations come from applying the 
total time derivative operator d/dt in cylindrical coordinates to the total mass and angular 
momentum associated with the annulus. In the case of angular momentum conservation, 
the derivative of the torque L follows from taking the limit of Aa —)■ 0. Combining these 
expressions we obtain 


aaS = 


da 




-1 


I dL I dL 

__ 

2ti da Tiafl da 


(46) 


where the expression on the right applies in the limit of Keplerian flow. Putting this expres¬ 
sion back into the continuity equation gives the evolution equation for the surface density: 


as _ _ 1 d _ _ 1 ^ I 1 

dt ada^ ada \ iraVL da 


(47) 


It remains to quantify the torque, either in terms of L or the corresponding radial drift a. 
Multiple sources of torque are just summed together, using both L (for viscosity) and a 
(for gravitational spreading). The l atter we inclu de in the form of Equation flTT|) . while the 
former we calculate next, following iPringlel (jl98ll ). 


The viscosity in an astrophysical ring relates the shear in the velocity flow {adQ/da, 
which is —3/214 for a Keplerian ring) to the torque on an annulus. 


L = ua ■ 27raS • a 


6a 


—STra^Hufl. 


(48) 


These expressions define the kinematical viscosity u. They may be easily related to the 
laboratory dehnition of dynamical viscosity by specifying a disk scale height and working 
with the force applied to the common surfaces between an annulus and its neighbors. 
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Putting the above result (Equation dH])) together with the expression for the radial 
drift rate from weak scattering (Equation yields the desired evolution equation: 


as 

'm 


i_a 

a da 


da ^ ^ 


(ctScigi-av) 


(49) 


In this expression, the first term on the right-hand-side gives the diffusive evolution of S 
from viscosity u. The second term corresponds to gravitational scattering. Under some 
circumstances, for example, a ring of monodisperse particles with a shallow surface density 
gradient and relative speeds higher than the escape velocity, we can derive an approximate 
viscosity for the gravitational interactions. 


TT ir 

‘^esc^-^K ^esc 


-^grav 


pr 


(v >vj. 


(50) 


This expression is comparable, to within a con stant of order unity, to similar expressions 
derived from the standpoint of kinetic theory (jShu fc StewartNl98511 . Since we derive the 
ring evolution directly from scattering theory, and do not seek a quantity of the form of a 
viscosity parameter, our theory differs slightly from the kinetic approach. In principle, an 
advantage of our approach is that it can accommodate migration when the distribution of 
ring particle mass covers a wide range. Below 11 13. 5. 2p we further discuss the connection 
between the evolution from the effective (collisional) viscosity and the contribution from 
gravity in our formalism. 


The solution of Equation (14^ requires boundary conditions. In general, the outer edge 
of the ring might expand into a vacuum (or be truncated artificially at the edge of a compu¬ 
tational domain). The inner boundary is more sensitive to the astrophysical setting. With 
a single central mass, material may be lost by accretion, thus S = 0 at its surface. For 
a circumbinary ring, the lack of stable orbits within a factor of at least twice the binary 
separation (for comparable-mass binaries) establishes the inner boundary of the disk. Pre¬ 
liminary simulations show that particles at the inner edge get pushed back into the ring, or 
are trapped in resonances there 11 14.11 below). We therefore provisionally set inner b oundary 
cond ition at the innermost stable circular orbit to be reflective, dS/da = 0 (see also [Pringle 


199l|). 


Together with specific values for z/ 11 13.3P and the drift rate from mutual weak scattering 
of ring particles 11 13.4p . these boundary conditions allow us to track the spreading of a ring 
from some initial state. 















20 


3.5.2. Connection between weak scattering and ring viscosity 


Ring viscosity and gravitational scattering have similar physical origins. In both cases, 
interactions between pairs of particles cause deflections which lead to a radial mass flux. 
There are also profound differences. Deflections in a viscous ring derive fror n physical colli¬ 
sions , which can stir particles up un til collisional damping becomes important fjCook fc Franklin 


1964J: iGoldreich fc Tremainelll978l ). Gravitational scattering can also stir particles, but only 


until they reach the escape speed. Beyond that speed, gravitational scattering becomes in¬ 
effective. Furthermore, gravitational interactions become irrelevant if the particles he inside 
the Roche limit. In both of these situations, the ring is completely driven by the viscosity 
from physical collisions. Otherwise, gravity and viscosity share in the evolution of the ring. 


To compare quantitatively how viscosity and gravity operate, we focus on the surface 
density flux terms (right hand side of Equation fl49lB . Under a set of simplifying assumptions, 
including constant viscosity u and surface density S in a ring of monodisperse particles, we 
write the fluxes in terms of the particle escape velocities Vgsc- 


4isc = 

Jgrav = aScigrav ~ GS vT^ ^2 ^ (^2^2 


(51) 

(52) 


where C ~ ~ 4. Thus, the particle velocity, as compared to the escape velocity, 

delineates the viscous and gravity regimes. Roughly, viscosity operates at n S> Uesc, while 
gravitational scattering drives the evolution at n <C Uesc- 

In a viscous disk, the effective viscosity formally depends on the random velocity with 
no restrictions except that v must small compared to the local Keplerian velocity. In the 
gravity dominated regime, the surface density flux has a maximum, 

Jgrav~GSVT^^ (^;<r;esc). (53) 


Physically this limit reflects our assumptions about gravitational interactions: particles with 
low relative speeds drifting by one another with impact parameters within are in the 
chaotic, orbit crossing regime, or possibly corotating. Here, we assume that these types of 
orbits do not lead to a net radial drift. In any event, if the speeds are below the escape 
velocity, then collisional damping also must be operating, keeping the ring dynamically cool 
on these small scales. 
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3.6. Illustrations of the theory 


To give examples of solutions to the evolution equation (jH]), we use a finite difference 
code that divides the surface density into 640 radial bins. (In runs described below, we 
doubled or quadrupled the number of bins to track the physics of any sharp edges that 
develop.) The initial surface density has a Gaussian profile in orbital distance, and a fixed 
total mass. In each run, we choose a particle radius and set the dispersion of the particles to 
be their escape velocity. These parameters, along with the surface density on the grid and 
its gradient (derived from finite differences) give the viscosity i^rad and cigrav needed to evolve 
E in a step-wise fashion. 

Figure [2] show numerical solutions for the evolution of a ring around a central mass. 
The curves in the Figure are snapshots of the evolution, corresponding to a viscous ring and 
to an inviscid ring that spreads from gravity. The two phenomena yield similar behavior in 
the ring at this delimiting velocity. 


For reference, we consider how the theory applies to Saturn’s rings. Inside the Roche 
limit, mutual gravity between small ring particles is irrelevant. Outside this limit, external 
influences (distant moons) stir the ring particles to speeds that exceed their escape velocity. 
For example, the A ring consists of small (centimet er- to meter-size) particles, with a scale 
height of roughly 10 meters fe.g.. lColwell et al.ll2009h . The orbital frequency is O ~ 10““^ s“^; 
typical random speeds, v ~ hfl, are ~0.1 cm/s. The escape velocity of a centimeter-size 
particle is roughly 0.001 cm/s; thus, viscosity dominates in the ring. 


In other planetary rings, conditions probably differ. Our focus turns to the Pluto- 
Charon binary where particle masses may have been large enough for gravity to have an 
impact on ring evolution. 


4. Circumbinary orbital dynamics 


Pluto and Charon, as partners in the solar system’s most prominent binary planet, 
pose an interesting challenge to our understanding of the dynamics of plan etary rings and 
satellite s. Circumbinary orbits differ from Kepleria n ones in important wavs flHeppenheimei 


1978 : Murray fc DermottI Il999l: Ihee fc Pealel 120061: Iheung fc heel l2013l: iBromlev fc Kenvon 
20151 ). In this section, we present an overview of these differences. 
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4.1. Orbital stability and resonances 

At distances comparable to the binary separation , ahin. circumbinary orbits are un stable. 
To help map out regions of instability around a binary, iHolman fc WiegertI fjl999l) numerically 
estimate a critical distance from the binary center of mass, Ocrit, that depends on the binary 
eccentricity (cbin) and mass ratio. Beyond Ocrit, most orbits are stable over many (10"^) 
dynamical times. Inside Ocrit, the gravity of the binary destabilizes orbits. Thus, every 
binary has an inner cavity that is cleared of ring particles. For Pluto-Charon in its present- 
day orbit (cbin ~ 0), this inner cavity lies inside of Ocrit ~ 2.0abin- With a ~ 2.2abin, the 
innermost moon Styx is just beyond this critical orbital distance. 


Instabil ities arise when reso n ances with the central binary overlap fjWisdomlll980l:lLecar et ah 


200l|). The IHolman &: WiegertI fll999l) expression for Ocrit gives only an approximate loca¬ 
tion of the outer edge of the central cavity around the binary. Some stable orbits exist 


i nside of Ocrit; there may be various unstable resona nces outside of it. Ponova &: Shevc 


f 2013h provide an excel lent illustration (seea also Wvatt et ah 1999 
Doolin fc Blundell 2011 ). 


Pichardo et ah 


lenko 


For the Pluto-Charon binary, with eccentricity near zero, resonance phenomena are 
not important except near the 3:1 resonance at ~ 2abin- When the binary has a modest 
eccentricity, Cbin = 0.1, however, unstable orbits form around the 5:1 commensurability. For 
other values of ebin, it is important to map the stability orbits f or each astrophysical setting. 


especially for prob lems which involve massive satellites (e.g. iPopova fc Shevchenkol 12013 
Cheng et ah 2014a). 


Figure [3] illustrates the effects of resonances. Initially, a massless satellite orbits outside 
of the 7:1 commensurability around an expanding circular binary. On its own, the satellite 
navigates the 7:1-4:1 commensurabilities as binary expansion causes them to sweep past. 
Near the 3:1 resonance (with a close to Ocrit), however, the satellite is ejected (Fig. |3l left 
panel). 

In contrast, a satellite that experiences sufficient collisional damping and diffusion sur¬ 
vives. In Fig. [3] (right panel), a satellite starts in the same orbit but is subjected to diffusion 
and random walks in a manner consistent with a monodisperse ring of 1-km particles stirred 
to their escape velocity in a rin g with surface densit y 40 g/ cm^ (see ^5.21 below). We mimic 


these effects in the manner of iBromlev &: KenvonI fl2013l. e.g.. Equation fl3TD and related 


discussion, above); damping comes from a slow and steady adjustment of the osculating 
Keplerian eccentricity, while diffusion is a sequence of random walks. 


With collisional damping and diffusion, satellites in the rings coexist with the destabi¬ 
lizing resonances. Strong damping (many small particles) allows particles to orbit as if the 
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instability did not exist; weak damping (fewer, larger particles) and higher binary eccen¬ 
tricity leaves satellites snsceptible to ejection. Intermediate levels of damping and diffnsion 
allow particles to become trapped in resonances and to move with them. In the example 
shown in the left panel of Fig. [31 stable satellites expanding in the 3:1 commensnrability 
reqnire damping times < 50 yr (more than 10^ satellite orbital periods). When all of the 
mass in the ring is in objects with radii smaller (larger) than a few km, the time scale for col- 
lisional damping (cf. Eqnation (ITT]) is smaller (larger) than this reqnirement. Tims, different 
mixtnres of particle sizes lead to different ontcomes for the evolntion of the Plnto-Charon 
binary and its satellites. 


4.2. Circumbinary orbits and most circular paths 


Ontside the resonances, circnmbinary orbits are stable and can be snrprisingly circn- 
lar. Despite the non-Keplerian gravitational potential of the binary, there is a set of orbits 
analogons to circnlar paths aronnd a single central mass. This family of ‘most circnlar ’ 


orbits consists of nested, circnmbinar y trajectories that do not cross flLee &: Peald 12006 


Youdin et al.ll2012l: iLenng fc Leell2013[) . Even when the orbit of the central binary is eccen¬ 


tric, t he most circnlar paths have eccentric motion ch aracterized by a ‘forced’ eccentricity, 
Cforce flHeppenheimerl Il978l: [Murray fc DermottI [19991) that preserves the nested and non¬ 
intersecting properties. 


Particles with random motion about a most circular orbit have a ‘free’ eccentricity, Cfree, 
which serves the same role as the orbital eccentricity for an orbit around a single central 
mass. Particles with free eccentricity collide with characteristic relative velocities CfreeW- 
Particles with eccentricity that damps to zero settle on most circular paths with a forced 
eccentricity. 


4 . 2 . 1 . Analytical theory of circumbinary orbits 


Most circular orbits appear as part of general circumbinary orbit solutio ns from ana¬ 
l ytica l the ory in the con t ext of the restricted three body problem. Following [Lee fc Peale 
fl2006[) and [Leung fc Leel (j2013|), we start with the potential 


$ = 


GM„ 


GM,. 




'^ + z^ + Rl 


2RRp cos A(j) + 2RR^ cos ’ 


(54) 


where Acf is the angle between the secondary and the massless satellite (ring particle) in a 
reference frame with the binary’s center of mass at the origin; in this frame, the satellite is 




























24 


at radial position R in the plane of the binary, and altitude 2 ; above it, while Rp and Rs 
are the orbital distances of the primary and secondary, with masses Mp and Mg. When the 
binary is circular (cbin = 0), and Rs are constants. 


The strategy of iLee &: Pealel (120061) is to expand the potential in terms of the angle 
cosines in Equation (l5T|) : for a circular binary, this expansion is 


$ ~ E $o/c cos(fcA0) 


[^bin 0, Z 0] 


(55) 


A:=0 


where flbin is the mean motion of the binary (flbin^ = G{Mp + Ms)/aur?)-i and our choice 
of time t hxes the orbital phase. Each coefficient rnay be expressed as a series in terms 
of R^/au^-n, where is the radius of a gu i ding c enter associated with a satellite’s orbit. 
Lee fc Pealel (120061) and iBromlev fc KenvonI (l2015l) list leading terms in these expansions. 


The next step is to express a satellite’s motion as excursions in radial, azimuthal and 
altitude coordinates {6R,6(j),z) about the guiding center at radius Rg. The result is a set 
of equations of motion that can be linearized in the excursion coordinates. When set up 
in the manner of driven harmonic oscillators, these equations can be solved with standard 
techniques. The natural frequencies Ke and ^ correspond to the free eccentricity and free 
inclination. The driving frequencies are the synodic frequency and its harmonics, along with 
the binary’s orbital frequency if the binary is eccentric. 

This prescript ion yields expressions for the satellite position in cylindrical coordinates 
( Lee fc Peale 2006 1: 


R{t) = 


^free '^e) ^force COs(l^gt) ~h ^ ^ Ck COs(/nCJsynt^ -|- 


k=l 




t H-^ sin(fi;e^ + i’e) + eforceSin(f2gf) + -—— sin(fca;syn^) + 




(56) 


(57) 


z{t) = iRg cos(z/it + 'ipi), 


(58) 


where f2g is the angular speed of the guiding center, cjsyn = f^bin —is the synodic frequency 
of the satellite relative to the binary; orbital elements are Cforce, Cfree and z, the forced 
eccentricity, the free eccentricity and the inclination, respectively, while the phase angles 'ipe 
and xjji are constants. The coefficients Ck and Dk are associated with the time-v arying part of 
the b inary’s potential. Terms not shown depend on the binary eccentricity (see heeling fc T^ee 
2 OI 3 I ). In the form given above, the solutions are accurate when eun^Rg/aun ^ 1- 
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To define the parameters in Eqnations fl56l) - fl58l) . we start by noting that the orbital 
freqnency of the guiding center comes from the time-average potential $oo (Equation (jSS])): 


= 


1 


00 


Rg dR 


GM 




1 + 


ii 

M 


3abi„2 45 (M3 + Mf) abin" 


4 R, 


+ 


64 M3 


Rp 


+ ... 


(59) 


the square root of which is the mean motion of the ring particle. Its epicyclic frequency is 


kI = R. 




dR 


- = 


Re: 


GM 


1 -A 

M 


3auJ , 135 (M3 + m 3) abin" 


4 Rp 


+ 


64 


M3 


Rp 


+ ... 


(60) 


A similar analysis based on the expansion of the potential in the ^(-direction yields the vertical 
excursion frequency, 


2 = 

z dz 


z=0,Rg 


GM 


1 + ^ 
M 


9abin" , 225 (M3 + m 3) abin" 
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Because and differ from the mean motion of the satellite, the free eccentricity and 
inclination have corresponding nodal precession. In the limit that the binary separation 
goes to zero, both Ke and Ui become the Keplerian mean motion. Next, the coefficients for 
the forcing terms at the synodic frequency are 
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in the notation of iLee fc Peald (120061 ). Finally, we have the forced eccentricity. 
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flHeppenheimerl Il978l: iMoriwaki fc Nakagawal 120041: iRafikovI 120131) . Iheung fc Led (12013l) in¬ 


clude this term {G^ in their notation), along with terms similar to Gk and Dk but reduced 
by a factor of the binary eccentricity. 


4-2.2. Examples 

Most circular paths, with no inclination or free eccentricity, describe the least dissipative 
trajectories for particles that suffer collisions. Small particles on these paths remain on 
them indefinitely without colliding. We construct these trajectories using the orbit solutions 
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in Equations fl56|l -fl58 |l . Because the fall-off with a^n/Rg is shallow for orbits that are 
modestly close to the binary, we also consider numerical search methods. We evolve many 


trial paths until we hnd the one that minimizes the radial excursion (e.g., lYondin et ah 


20121) ■ This procedure is slow, as it requires numerical orbit integration with each trial. For 


the integrations, we use a 6^^-order symplectic time step algorithm in the N-body component 
of our Orchestra code (IBromlev fc Kenvonll2006l. l2011a|) . 


Figure m provides an illustration of a particle on a most circular orbit around a binary 
with Cbin = 0. It shows a direct comparison between the analytical theory and a most circular 
orbit found by numerically searching for the trajectory with minimum radial excursions about 
its guiding center. The small radial excursions of a particle on a most circular orbit have a 
“W” shape prohle over a synodic period, showing (i) maxima when the satellite is aligned 
with the binary and closest to the secondary, (ii) minima when the satellite is perpendicular 
to the binary axis, and (iii) local maxima when the satellite is closest to the primary. 


In the circular binary’s frame of reference, the orbit of a satellite is a closed, egg-shaped 
loop. All such hgures are nested. The further away from the binary, the more circular the 
hgure. Particles on most circular orbits traverse these paths on their own synodic period 
with respect to the binary. While their speeds may be different, their phases align in a 
way th at all neighbori n g part icles get coherently jostled by the binary. None of their paths 


cross. 


Pichardo et al.l (120051) provide an excellent description of these orbits in terms of 


‘invariant loops’, closed cur ves that vary in response to the binary’s potential (see also 


Georgakarakos fc Eggll 120151) . 


If the binary has some eccentricity, most circular orbits are eccentric. Satellites following 
these paths have some epicyclic motion about the guiding center, with eccentricity set by Cforce 
in Equation (1641) . Figure [5] provides an illustration, using a small binary eccentricity to allow 
the radial excursion from the forced eccentricity to roughly match the radial displacements 
oscillating at the synodic frequency. With zero free eccentricity, the radial excursions (~ 
.RgCforce) are steady, with no overall change in their maxima over time. 

Figure [5] also shows the effect of free eccentricity. Although the free eccentricity has 
nodal precession, the forced eccentricity does not. Thus, the radial excursions relative to 
a purely circular orbit about the binary center-of-mass vary continuously, with a period 
27r/|f2g — Ke\- Adopting a free eccentricity equal to the forced eccentricity for the example 
in Figure O the total eccentricity cycles between zero and 2eforce- Thus orbits with no free 
eccentricity (cfree = 0) have the smallest maximum radial excursion. More importantly, these 
orbits do not intersect, offering a place in phase-space for many small particles to coexist 
without colliding with each other. 
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4.3. Summary 


We conclude this section with its main implication for ring dynamics. The family of 
most circular orbits represents the dynamically cold state of a ring orbiting a central binary. 
Collisional damping causes particles to settle onto one of these orbits, just as for circular 
orbits around a single central mass. As long as the ring is not globally self-gravitating, the 
most circular orbits represent reference frames in which the local evolution of the ring takes 
place. Signihcantly, these orbit solutions exchange no net torque with the binary. Thus, the 
theory for ring dynamics developed for a purely Keplerian potential applies to circumbinary 
orbits, as long as the guiding centers are at distances where long-term stable orbits are 
possible. 

Here we take advantage of these results by adopting reference frames that track most 
circular paths. Then we may perform calculations such as solving for the radial evolution 
of a circumbinary ring (Equation simply by setting the orbital distance variable a to 
the guiding center radius. There are still differences, although they are typically small. For 
example, the viscosity in Equation (fT7|) depends on the orbital period, which is reduced by 
a factor of 3/8 x {^/M){aun/Rg) in units of the binary separation. At the orbit of Styx, this 
correction is about 1%. 


The strategy of working in most circular frames fails at distances close to binary or 
at the locations of resonant instabilities. Howe ver, the effects of destabilizing multiple 
resonances (e.g., IWisdomI Il980l: iLecar et al.l l200l[l are avoided if particles are collisionally 
damped. Furthermo re, as in Figure 3l da mping and diffusion may facilitate resonant trap¬ 
ping and migration fjWard fc Canupll2006|) . processes that can enable a ring to expand as its 
central binary evolves. Preliminary estimates suggest that the damping time expected for a 
sea of 1 km particles — roughly 10 yr — is more than sufficient to mitigate the effects of the 
instabilities outside of the 3:1 commensurability with the binary. On their own, signihcantly 
larger objects, with radii exceeding 5 km or so, may not be able to orbitally damp on their 
own, but instead require a population of smaller bodies to keep them stable. 


5. The Pluto-Charon system 

In this section, we apply the theory of circumbinary ring dynamics and evolution to 
Pluto, Charon, and a ring where satellites form. Aside from deriving quantitative estimates 
for time scales and rates, we numerically solve for the evolution of possible ring conhgurations. 
We begin with a consideration of the binary itself, since its evolution certainly affects the 
ring. 
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5.1. Formation and evolution of the binary 

1153 km and Rc. = 603.6 km and masses Mp 


Pluto and Charon have radii -Rp 
1.305 X 10^^ g and Me. ~ 1.520 x 10^^ g fiBuie et al.ll2006l: IPerson et al.ll2006l: iTholen et ah 
2 OI 2 I: iBrozovic et al.N2015[) . Their orbit h as separation ap(^ 16.97 i?p, eccentricity epc 


consistent with zero, and is tidally locked flBuie et al.ll2012l) . The binary orbit relative to 
the orbital plane of the major planets is significantly inclined, with ipc ~ 96°. 

An impact origin for the Pluto-Charon is plausible, both in terms of the collision out¬ 
comes as compared to the b inary’s present configuration and the likelihood of such a collision 
in th e outer solar s ystem ( McKinnonI Il989l: ISternI Il992l: ICanupl 120051: iKenvon fc Bromlev 
2014). Simulations f Canup 20051) suggest that as much as 3 x 10^^ g — about 0.02% of 
the binary’s mass — ends up as debris around the binary. This material initially lies on 
yery eccentric orbits extending from distances of 5-10 R-p to 30 Rp. If the orbits circu- 
l arize, the debris probabl y forms a compact ring within 15-30 Rp of the center of mass 


( Kenyon &: Broml^ 2014 ). 


Kenyon fc Bromley! (120141) propose that collisional evolution of this ring leads to the 


formation of the four known satellites and perhaps other, smaller satellites at larger distances 
from Pluto-Charon. The known satellites — Styx (at an orbital distance of approximately 
37 Rp), Nix (43 Rp), Kerberos (5( 3 Rp) and Hydra (57 Rp) — have a combined mass in the 
range of 3 x 10^^ g to 3 x 10 ^^ g ( Buie et ah 2006 : Tholen et ah 2012 : Brozovic et ah 2015 


Showalter fc Hamilton! 120151) . The combined mass is consistent with formation from debris 
left over from a collision which produces the Pluto-Charon binary. The orbital distances 
require expansion of the ring. 


In an impact scenario, Pluto and Charon haye an initial orbital separation of 

flpc(^ = 0) ~ 5i?p 


(65) 


(jCanupll2005f) . The separation grew in time through exchange o f angular momentum that 
depends on the state of matter within the two partners (e.g., iBarr fc Colhn.sl 120151) : an 
idealized theoretical rate describing the expansion is 


da 


PC -11/2 

= KOr 


dt 


"'PC 


where 


K = 3- 


k Mr 


G{Mp + Mc) ] 
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( 66 ) 


(67) 


QMp [ 

with Love number k = 0.058 and a constant tidal dissipation function Q = 100 (e.g., 
Farinella et ah 1979 : Dobrovolskis et ah 1997 : Peale 1999; Cheng et ah 2014a ). The solution 
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to these equations thus describes the orbital separation as a function of time: 

apc{t) = niin I [l3Kt/2 + apc(O)^^^^] ^ ^ , 17i?p|. (68) 

In this model, the binary expands by a factor of two within 10® yr and reaches its present-day 
separation within 3 Myr, when Pluto and Charon become tidally locked. 


The initial eccentricity affects the tidal evolution, typically by speeding up the expan¬ 
sion. However, for modest eccentricity, the effect is weak, and in any event, the eccentricity 
damps relatively quickly: 


1 depc{t) 


epc 


dt 


F - 13/2 
- —KOpc 


(69) 


where the coefficient F is 0(10) fjCheng et ahl l2014a|) . The damping time for an initial 
eccentricity of order unity in this case is of the order of 1000 yr. Hence, we assume that the 
binary eccentricity epc is zero unless otherwise specihed. 


5.2. The Pluto-Charon ring: early conditions 


To describe the post-impact conditions in the Pluto-Charon system, we dehne a hducial 
model with a binary separation of Obm = 5 i?p, and circumbinary ring radius of Oring = 20 Rp. 
We assume that the mass of the ring is 


Mring = .Yri„g ■ 3 X 10“ g 


(70) 


where Wing > 1 is a rough bou nd for our models in which the ring must contain enough mass 


to make the known satellites ( Kenvon &: BromlevI 12014) . We assume that this material has 
a Gaussian surface density prohle with standard deviation of Aa = 2 Rp: 


E(f = 0) = Scexp [-(a - /2Aa^]^ 
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(72) 


where Eg- is a constant relevant only to this Gaussian model. 


To estimate plausible characteristic scales, we consider a monodisperse population of 
ring particles with r = 1 km and p = 1 g/cm®, and a characteristic surface density of 
40 g/cm^. The number of particles in the ring is 
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7 X lO^X 


ring 


1 km 
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(73) 
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The speeds of individual particles that enter into ring dynamics are 

1/2 


= 0.21 


20i?p 


= 16 


Vpec, — 75 


- 1/2 


r 


20i?p 


1 km 
r 

1 km 


km/s 

[orbital speed] 

(74) 

cm/s 

[escape speed] 

(75) 

cm/s 

[Hill speed]. 

(76) 


The latter two speeds are characteristic of gravitational stirring in the ring. 

An important con sideration for a compact ring is whether it is stable against its own 


gravity flToomrel 11964) . To avoid gravitational clumping of particles within the ring, the 


particles must maintain a minimum random velocity in addition to their orbital motion. For 
stability, the characteristic minimum speed is 


vq = 0.94 


a 

3/2 

E 

20 Rp 


40 g/cm^ 


cm/s. [minimum stability speed]. 


(77) 


In a ring of identical icy particles, this condition translates to a radius of about 10 meters if 
they stir each other to their escape speed. Since the ring expands to larger radii and lower 
surface densities, the minimum speed and radius both fall as the ring evolves. Furthermore, 
if the impact debris contains larger embedded fragments, it will settle to higher speeds, 
comparable to the escape velocity of those fragments. Thus the bulk of the ring material 
can be composed of smaller particles and yet the ring will remain gravitationally stable. 

In addition to the binary separation and the orbital radius of the ring, there are several 
other useful length scales: 
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(78) 

(79) 

(80) 


Finally, results from §2 provide relevant time scales. Dynamical times are 
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d [binary orbital period] 
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Ring evolution time scales are 
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(83) 

(84) 

(85) 


where the migration time Tmigrate ~ o/hgrav, and Tgpread is (a/2)^/i/, giving a rough time 
scale for the ring to double its size by viscous diffusion. The gravitational spreading time 
scale is a crude estimate in the limit of a dynamically cold disk; dynamically hot rings, with 
scale heights that greatly exceed the particle Hill radii, will experience slower evolution by 
gravitational scattering. The numerical value for the viscous diffusion time carries the as¬ 
sumption that the random component of the particle velocities (governing the ring viscosity) 
are given by the escape velocity from the particle surface. The surface density E refers to 
the local surface density in the ring. 


For a ring of monodisperse particles with random velocities comparable to particle escape 
velocities, damping is generally much faster than spreading. For small particles (r < 1 km) 
in low-mass rings (E < 40 g/cm^), evolution time scales are comparable to the expansion 
time of the binary (10^-10® yr). Large particles or large ring mass can shorten the time 
for the ring to spread. However, it is hard to reduce the spreading time signihcantly unless 
particle sizes rea ch 50-100 km. The small nurnber of ring particles then makes spreading 
unlikely (see also lPeale et al.ll201ll: ICheng et al.ll2014m . 


If the ring consists of particles with a broad range of sizes, viscous diffusion can drive 
a rapid expansion of the ring. If a few large particles (radius riarge) are embedded in a sea 
of small ones (radius r smai i), gravitational stirring by the large particles drives the random 
velocities of the smaller particles to the escape velocity of the large particles. If the large 
particles have sufficien tly large escape velocities, the formal spreading time is very short 
I Goldreich et ah 2004): 
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( 86 ) 


This expression assumes that the rate of collisional damping for the small particles 
is small compared to the stirring rate of the large particles. For a swarm of 1 m par¬ 
ticles, the likely damping time of 0.01 yr (Equation is smaller than this spreading 

time. However, rapid damping promotes the growth of small objects into larger ones (e.g.. 
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Kenvon &: Bromlevll2014j) . slowing the rate of damping and accelerating the rate of spread¬ 
ing. Becanse the damping time grows mnch faster with a than the spreading time, it seems 
plausible that a swarm of particles with sizes ranging from 1 m to 1-10 km could spread 
radially on time scales ranging from a century to many millenia. 

While a two-population model for the ring is simplistic, it captures the relevant physics 


Wetherill & Stewart 

1989; 

Ida & Makino 

1993; 

Goldreich et ah 

2004) 


the early evolution of the debris is likely dominated by destructive collisions of particles 
on very eccentric orbits. For typical collision velocities of 0.1-0.2 km/s (Equation (El])), 
center-of-mass collision energies of 10^ — 10^ erg/g destroy particles with r < 5-10 km 
flBenz. fc Asphaud Il999l: iLeinhardt fc StewartI 120121: iKenvon Sz Bromlevi 120141) . Large self¬ 
gravity saves larger particles from signihcant mass loss. Rapi d collisional damping likel y 
saves much smaller particles, r < 1-10 m, from destruction (IKenvon fc Bromlevi 120141) . 
Fragmentation likely converts intermediate-mass particles, 1-10 m < r < 5-10 km, into 
smaller particles which damp rapidly. Thus, two populations of particles — small particles 
with most of the mass and a few large particles with signihcant self-gravity — are plausible 
components of a circularized ring of particles around Pluto-Charon. 


5.3. Evolution of the Pluto-Charon ring 


To investigate the evolution of monodisperse and two-population ensembles of particles 
in a ring orbiting Pluto-Charon, we solve the radial diffusion equation (Equation dH])) using 
a hnite difference code (see also iBromlev fc Kenvonll2011aj) . The code tracks the surface 
density S in 0(10^) radial bins spanning a = 5 Rp out to 120 Rp. In each bin, we calculate 
the viscosity (Equation ffT7|) i and the radial mass hux from gravitational scattering (Equa¬ 
tion fH2|) ). For each time step, the hnite diherence code updates t he surface density from 


the v i scosity and radial rnass hu x. Tests on standard problems (e.g., Lvnden-Bell fc Pringlel 


19741: iBromlev &: Kenvonll2011a|) verify the accuracy of our solutions to better than 1% over 
10® - 10^ yr. 


In applying our code to the circumbinary environment, we adopt the results of ^ 
Damping and stirring of particles, along with radial transport, take place relative to most 
circular orbits. Hence, our use of orbital distance a is synonymous with most circular orbit 
positions characterized by the orbital radius of the guiding center Rg. We also assume that 
resonant excitations and instabilities within the ring do not contribute signihcantly to the 
evolution. In our hducial ring with 1 km particles, the damping time is rapid enough to 
suppress resonant-driven excitations. These idealizations still allow us to reach our goal of 
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understanding how the ring might spread as the central binary expands. 


To establish the range of possibilities for the expansion of a compact ring around a 
young Pluto-Charon, we choose from a palette of models to test the effects of particle size 
(r), velocity (n), and total mass in the ring (Mring). Our initial conditions assume a Gaussian 
surface density prohle (Equation flTTD l with a = 20 Rp, Aa = 2Rp. We start the binary at 
an orbital separation of 5 Rp, and have it expand on a time-scale of ~ 10® yr, according to 
Equation 


An additional model parameter is the choice of the inner boundary condition. We can 
set the inner boundary to be “flow-through” under the assumption that all m aterial reachin g 
the inner edge is lost, presumably to ejection or accretion by the binary (e.g., |Pringlelll99lh . 


However, in most of our trials, the inner boundary is reflective at the edge of the stability 
zone. We enforce this boundary by setting up an artihcial outflow that is constant inside 
the ring edge and falls off exponentially over a scale length of 2.5% of the orbital distance of 
the edge. This algorithm prevents a strong discontinuity from building up at the inner edge 
of the ring. Even with this approach, the gradient at the edge can be large and can lead 
to numerical instability. Thus, we add artihcial viscosity in the vicinity of the edge in cases 
where gravitational scattering strongly dominates over viscous spreading. 


For all models we set S to zero at the outer boundary, although in most cases the ring 
never reaches the outer edge of the computational domain. 


The following list summarizes our models. We use the approximate time for the ring 
to reach 60 Rp, just beyond the present-day orbit of Hydra, as a benchmark. If this 
time is within ^ 10^ yr, then the ring may spread before large satellites grow within it 
flKenvon fc BromlevI 120141 1. Otherwise, some evolution of particle size — or some other 
means of migration — is required as the ring expands. 


• Baseline model. Our starting point is a ring with total mass 3 x 10^° g (Wing = 1) 
made up of kilometer-size particles stirred to their escape speed (~75 cm/s). Figure [6] 
illustrates the spreading of this ring. This ring spreads slowly compared to satellite 
growth times, extending past 60 Rp in a few times 10® yr. 

• Massive rings. If the ring is more massive than the baseline model, the spreading 

time is reduced (Tgpread 1/S 1/Wring). Figure [7] gives an example for a ring with 

Mring = 3 X 10^^ g (Wring = 10). In this model, the ring reaches 40 Rp in a 10^ yr and is 
beyond 70 Rp by 10® yr. Compared to the estimated growth time scale of 10^ — 10^ yr, 
this expansion rate is still slow. 

• Massive rings with loss through the inner boundary. In previous models, we use re- 
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flective inner boundaries, since material at the inner edge is presumably pushed back 
into the ring by the binary. Here we explore the possibility that mass instead flows 
through the inner edge and is ejected or accreted by the binary. Figure |H] shows an 
example. It shows the baseline model, as well as a model with a ring that is ten times 
more massive, both with and without loss of mass through the inner boundary. The 
effect of the mass leakage is a loss of about 20% of the total mass, a flattening of the 
surface density profile, and a slightly less extended disk. 

• Growth of particles. We start with the baseline model, but with particles of radius 
= 10 m. These particles grow monotonically to a maximum radius of rmax = 3 km 
on a time scale of rgrow according to 

r{t) = rmin + (rmax “ rmin) [l “ exp(-f/rgrow)] • (87) 


The random velocity is set to be the instantaneous escape speed of the particles. 
Figure [9] shows the result. It differs from the baseline model because the spreading 
starts off more slowly and becomes more rapid only at the late stages, as r grows. Even 
so, the ring extends to only 50 Rp in 10^ yr. 


Stirr ing by embedded satellites. Guided by the physics of planet formation flKenvon fc Bromlev 


20041 ). we Anally consider a scenario in which a population of small (r = 1 m) particles 


is stirred by a few larger ones (r = 10 km). We include the effect of these larger bodies 
by stirring the smaller ones to a random speed of half of the Hill velocity of 10 km ob¬ 
jects (fixed at its value at a = 20 Rp). This choice acknowledges that the larger bodies 
are rare, and that damping among the smaller bodies prevents them from being stirred 
to the escape speed of the embedded satellites. The result is presented in Figure [TOl 
The ring expands to Hydra’s present-day orbit in just over 1000 yr, and by 10“^ yr, it 
has spread to well over 100 Rp. 


To summarize our results, monodisperse populations of ring particles stirred to their 
own escape velocity, have a spreading time that decreases with particle radius and mass of 
the ring (Equation flS^ i. For particles with radii larger than a few kilometers and for ring 
masses larger than ~ 10^^ g, the ring spreads to Hydra’s orbit within 10® yr. Allowing for 
particle size to grow from some small value (r ~ 10 m) only delays the spreading. 

In models inspired by planet-forming disks, with small particles and embedded satellites, 
the ring can spread quickly (Equation fl5B]) j. Even a low mass ring (3 x 10^° g) can spread 
out past Hydra on a time scale of ~ 1000 yr. A weakness of this scenario is that the stirring 
is kept at a constant level throughout the ring as it spreads. For this scenario to be possible, 
embedded satellites must either migrate within the disk, or form in the outer regions of the 




35 


disk as it spreads. This latter possibility seems plausible, si nce formation times can b e short 
if the ring should dynamically cool by collisional damping flKenvon &: Bromlevll2014j) . 


6. Conclusions 

We consider the evolution of a ring of particles around a central binary, with a focus on 
the Pluto-Charon system. Compared to earlier theo ries for evolving disks of particles (e.g.. 
Cook fc FranklinI 11964 iGoldreich &: Tremainelll978l) . we add two new features: (i) gravita¬ 
tional scattering, which we quantify in terms of weak scattering theory, and (ii) the impact 
of the circumbinary environment on ring evolution. Gravitational scattering by massive 
particles can contribute signihcantly to the expansion of the ring. As for the circumbinary 
enviro nment, aside from defining an orbital plane and creating an inner edge by destabilizing 
orbits flHolman fc WiegertI 119991 1. the impact of the binary is surprisingly small. 


Our main results are listed here, followed by a summary of our conclusions. 


Ring evolution. Our analysis of weak gravitational scattering provides a natural 
extension to the viscou s ev olution of particle rings and disks. The expressions of 
Cook fc FranklinI (119641 ) and IGoldreich fc Tremaine! (119781 ) hold until gravitational in¬ 
teractions start to become important. The demarcation between the viscous and grav¬ 
ity regimes is the ratio of the random speeds of particles to their typical escape speed. 
If the random speeds are large relative to Vesc, the particles do not feel their mutual 
gravity — evolution is viscous. When particles orbit beyond the Roche limit and 
have random speeds comparable to or smaller than their escape velocities, they stir 
one another gravitationally. Our results allow a smooth transition between these two 
regimes. 

Most circular orbits. The circumbinary environment r egulates the inner edge of a 


Wisdomlll980l: Murray &: Dermott 


ring and can destabilize orb i ts near close-in resonances (e .g., 

19991: iPichardo et al.ll2008l: iPopova fc Shevchenkol 120131) . Although the binary’s time 
varying gravitational potential modihes orbits in the rest of the ring, the family of 
most circular orbits provides reference frames similar to t he circular orbits around a 
single central mass (iLee fc Pealel 120061: iLeung fc heel 120131) . Dissipative disks or rings 
settle on to these orbits, and perturbations (e.g., from gravitational stirring) cause 
particles to acquire free eccentricity and inclination relative to these frames. Thus, we 
can calculate r ing evolution in the same way around a binary as we do around a single 
planet or star (iBromlev fc Kenvonll2015h . 
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Secular excitations. In initiating the research reported here, our exp ectation was 
that s tirring by the binary might help the viscous spreading of the ring (e.g., 


Moriwaki &: Nakagawa 


2004 iPaardekooper et al.l 120121: lRafikovll2013l: iLines et ahll2014j) . However, collisional 
damping causes particles to settle on most circular paths, even in the presence of a mas¬ 
sive axisymmetric disk and distant massive perturbers. Thus, the central binary is not 
necessarily a source of stirring or torque exchange. We explore the implications for the 
early stages of circumbinary planet formation in a separate paper flBromlev fc Kenvon 
2015h . 


Resonance effects and instabilities. The binary can still make a mess of orbits. 
In units of the binary separation, the Pluto-Charon sat ellites are closer in than known 


circumbinary planets from NASA’s Kepler mission fe.g. iDovle et al.ll201ll:IWelsh et ah 


2 OI 2 I: lOrosz et ahl l2012alJbl) . And, unlike most of the exoplanets, the Pluto-Charon 
moons orbit inside of suspected resonance-driven instabilities. If the moons emerge 
early in a compact ring, then they can get caught up i n res onant orbits that grow 
along with the tidally evolving binary fjWard fc Canupl l2006l) . While they may be 
transported in the resonances, they may just as easily get ejected, as our own sim¬ 
ple experiments demonstrate (Figure [3]). For the Pluto-Charon system, coordinating 
resonance transport in a way that preserves all of the moons does not seems possible 
f Cheng et ah 2014a '). 


Collisional damping can stabilize resonant orbits. Our simple simulations also 
demonstrate how collisional damping and diffusion can protect satellite orbits fro m ejec¬ 
tion (see also Ward Sz Canup 2006 : Lithwick fc Wu 2008 : Cheng et ah 2014a '). Thus 
a ring of small particles can survive the sweep of resonances during tidal expansion. 
Similarly, as long as a substantial reservoir of small particles is in place, we speculate 
that between collisional damping among the small objects and dynamical friction be¬ 
tween these small bodies and the large ones, the entire population is protected agains 
unstable resonances. 

Evidence for resonant migration. Our simple simulations also d emonstrate how 
collis ional damping and diffusion can facilitate resonant migration fjWard fc Canup 
20061) . Small particles become trapped near resonances and then get swept up in them 
as the binary expands. Figure [3] illustrates the phenomenon and highlights the role of 
binary eccentricity in the process. If a reservoir of small particles is maintained during 
the binary’s tidal evolution, small moons might also participate in this migration. It is 
beyond our scope here to explore these possibilities in detail, but our anecdotal results 
are promising. We plan more detailed simulations to test this possibility. 
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On the basis of these ideas and our numerical simulations, we infer the following con¬ 
clusions for the Pluto-Charon system: 


Starting with a compact ring around Pluto-Charon with mass comparable to the ob¬ 
served moons, viscosity and gravity can spread the ring to the moons’ present loca¬ 
tions in less than 10^-10® yr. These time scales are comparable to the tidal evolution 
of the binary’s semimajor axis, but are substantially longer than the growth time of 
10-100 km satellites (~ 10^ yr). 

The spreading of the ring can be comparable with, or faster than, the growth time if 
(i) the ring has small (1-meter) particles stirred by large (1-10 km) bodies, or (ii) the 
mass of the ring is large (> 10^^ g). In scenario (i) large objects can probably form as 
the ring spreads. In both scenarios, the ring spreads quickly to distances well beyond 
Hydra, the outermost known moon. 


Reson ance-driven instabilities remove isolated satellites or moons (e.g., iCheng et ah 
2014al) . This phenomenon is complicated by the tidal evolution of Pluto-Charon, 
which causes commensurabilities to sweep through the system. However, collisional 
damping and diffusion can make satellites immune to resonance-driven instabilities 
and can also facilitate trapping in a resonance. In the latter case, ring material can be 
tr ansported outwa r d alon g with the resonances, as envisioned for full-fledged moons 
by Ward fc Canup ( 2006 1. 


Together, the spreading mechanisms and interactions with resonances suggest poten- 
tial growth modes for th e satellites. In situ formation can occur quickly (10^ yr; 


Kenvon fc Bromlevi 120141 : the time needed for spreading the ring can be comparably 
fast or slow (10® yr), leading to these possibilities: 

i. The spreading might occur before the expansion of the binary. Then the sweep of 
resonances might eject or collisionally destroy the largest satellites. They would 
have to re-emerge out of the debris once the sweep was complete. This scenario 
would require a ma ssive initial ring, as particles could be lost in the process. 
flKenvon et al.l I2014J describe a similar picture for Fomalhaut b.) Simulations 
suggest that there was enough mass in the ring initial ly to sustain a ma ss loss 
of over 99% and still have enough to make the moons f Canup 20051. 2011 1. The 
result could be an extended disk and small as-yet-undetected moons beyond Hydra 
(IKenvon fc Bromlevi 12014 1. 


ii. If the ring spreads slowly, then particles in it could get caught up in one or more 
resonances with the expanding binary. In situ formation during the expansion 
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is possible, but a reservoir of small particles (perhaps replenished by collisions) 
would be essential to prevent the instabilities from leading to ejection. Preliminary 
numerical experiments (Figure [3]) show that material might be trapped in several 
different commensurabilities, depending on the binary eccentricity and the particle 
damping/diffusion. 


The next step is to incorp orate the physics of ring spreading into Orchestra, our hy¬ 
brid N-body-coagulation code (iBromlev fc KenvoiJ 1201 fall . Then, we can track the detailed 
interplay between ring spreading and satellite growth. The N-body part will be essential 


to sirn ulate resonant-driven instabilities that can clear out orbits in the ring (ICheng et ah 


2014b|). How these instabilities play out will depend on a number of factors including the 
effectiveness of orbital damping, ar id the time-varying eccen tricity of the binary, since cir¬ 
cularization brings some stability fiHolman fc WiegertI Il999l) . Our code will enable us to 
determine how collisional damping, gravitational stirring, resonance-driven instabilities, ring 
spreading and tidal evolution can work in concert to make the binary’s intricate system of 


moons. 


Meanwhile, our models suggest that the expansion of a narrow ring around the Pluto- 
Charon binary is plausible. The extent of the expansion depends on the details. The presence 
of small moons beyond Hydra’s orbits would suggest a more massive disk with small particles 
stirred by larger ones. Then, we expect 3-5 satellites with radii of a kilometer or so, on 
circular orbits beyond 60 i?p ( Kenvon &: BromlevI 2014 ). Ring evolution models with less 
mass and less stirring leave the ring truncated at Hydra’s orbit. Even so, scattering by 
the growing moons could place some satellites beyond Hydra, although their eccentricity is 
expected to be signihcant. New Horizons, with its unprecedented view of the Pluto-Charon 
satellite system, will test these scenarios. 


We are grateful for the comments of an anonymous referee that resulted in a number of 
improvements to our presentation. We thank NASA for support through the Astrophysics 
Theory and Origins of Solar Systems programs (grant NNX10AF35G) and through the 
Outer Planets Program (grant NNX11AM37G). We are also grateful to NASA for a generous 
allotment of supercomputer time on the ’discover’ cluster. 
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orbital distance (Rp) 


Fig. 1.— A schematic illustrating the layout and evolution of the Pluto-Charon satellite 
system. The binary is on the left. The satellites appear at their orbital distance from the 
center of mass, with symbol sizes scaled according to their approximate radii. The region 
inside of 3Rp is the Roche zone. Within this region, particles are held together by their own 
material strength; collisions dominate the particle dynamics. The light gray shaded region 
designates the inner cavity, where most planar orbits are unstable. The gray curves extending 
to large orbital distances locate commensurabilities (as labeled) as the binary expands. The 
brown dots at the lower left indicate the ring of debris after the binary formed, while the 
other clusters of brown dots refer to fast and slow modes of ring spreading. 
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Fig. 2.— A comparison between the effects of viscosity and gravitational scattering in the 
spreading of rings. The astrophysical context is Pluto-Charon, at some early time when the 
binary is compact (orbital separation opc = 5 Rp) and surrounded by a 3 x 10^° g ring. Each 
scenario begins with the same surface density prohle in the ring, but differ in the size of the 
ring particles, as indicated. In each case, the particles are stirred to their escape velocity. 
The dark curves spread from gravity only; light curves spread solely from viscosity. At this 
choice of particle speed, the two mechanisms should be similar. At higher speeds viscosity 
will dominate; at lower speeds gravity will drive the spreading. 
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Fig. 3.— The orbital evolution of a satellite around an expanding binary. The left panel 
shows a circular “Pluto-Charon” binary expanding from an orbital separation of 5 R-p (lower 
black curve). The dashed lines indicate the location of the 3:1, 5:1, and 7:1 commensurabili- 
ties expanding with the binary. A satellite on its own (black dots on the far left) orbits at a 
radial distance of 20 Rp until it is overtaken by an instability near the 3:1 commensurability 
at about 400 yr. Then it is ejected from the system over a period of hundreds of orbits. A 
satellite that experiences collisional damping and radial diffusion (gray dots) survives the 
instability and “surfs” the 3:1 commensurability as it expands. The damping and diffusion 
parameters are drawn from hducial values in ^5.21 (r = 1 km, v = Vgsc, F = 40 g/cm^). The 
right panel is similar, but with an eccentric binary (cbm = 0.1). In this case, an isolated 
satellite (black dots) is unable to navigate the 5:1 commensurability. A damped, diffusing 
satellite survives and rides resonances near the 4:1 commensurability. 
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Fig. 4.— The radial excursion of a ring particle on a most circular circumbinary orbit for a 
binary with Cbm = 0. The darker curve is from points calculated with a simulation starting 
from a best-£t to a minimum excursion orbit, while the lighter curve is from analytical theory 
(the two curves follow one another very closely). The binary separation and orbital distance 
correspond to the satellite Nix around Pluto-Charon. 
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Fig. 5.— The radial excursion of ring particles in orbit about a binary with Cbm = 0.01. The 
upper dark curve is for a most circular path, while the lower dark curve (displaced vertically 
for clarity) is for a particle that has equal measures of free and forced eccentricity. In time, 
the growth of the lower curve will continue; it will produce a beat patte rn with a maximum 
ampl itude equal to twice that of the most circular curve (e.g.. Fig. 3 of iBromley fc Kenyon 
20151) . The light shaded curves are from analytical theory. The orbital elements are chosen 
to be similar to Kerberos around Pluto-Charon. 
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Fig. 6.— Spreading of a circumbinary ring. The curves show snapshots of the surface density 
as a function of orbital distance from Pluto-Charon. The rings spread from viscosity and 
gravitational scattering as the binary expands. The binary pushes the inner edge of ring 
outward, keeping it at a distance of twice the binary separation. 
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Fig. 7.— Spreading of circumbinary rings of two different masses as in the previous figure 
but with prohles shown only at 10^ yr. The rings initially differ only in their surface density; 
the outcomes show how the increased surface density drives the spreading. 
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Fig. 8.— Spreading of a massive circumbinary ring with mass loss at the inner edge. The 
curves show the effect of mass loss at the inner edge of the ring, which expands along with 
the binary. 
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Fig. 9.— The spreading of a circumbinary ring with growing particles as in Fignre |6] except 
the mass grows (along with the particle speeds, set by the escape velocity) from 10 m to 
3 km on a time scale of 10^ yr. The growth allows for modest spreading beyond the baseline 
model with particle radii fixed at 1 km. 
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Fig. 10.— Rapid spreading of a circumbinary ring of small particles with embedded satellites. 
A population of meter-size objects stirred to 0.5 times the Hill velocity of 10-km objects 
(evaluated at 20 Rp) spreads rapidly, as seen in the sequence of curves. In this case, material 
can flow through the inner edge of the ring. While this model is an idealization, it shows 
how effective viscous spreading can be when the relative velocities are pumped up. 




